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Abstract 

(^ , Critical transitions are sudden, often irreversible, changes that can occur in a large variety of complex 

Cn ' systems. Signals that warn of critical transitions are highly desirable, but their construction can be 

impeded by limited availability of data. We propose a method that can significantly reduce the amount 
of time series data required for a robust early warning signal by using other information about the 
system. This information is integrated through the framework of a generalized model. We demonstrate 
the applicability of the proposed approach through several examples, including a previously published 
fisheries model. 
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Fisheries, coral reefs, productive farmland, planetary climate, neural activity in the brain, and finan- 
cial markets are all complex systems that can be susceptible to sudden changes leading to drastic re- 
Q^l organization or collapse. A variety of signals based on analysis of time-series data have been proposed 

that would provide warning of these so-called critical transitions. The amount of data these approaches 
require, however, can be prohibitive, especially in ecological contexts. We propose a new method for 
^ ' calculating early warning signals that can significantly reduce the amount of data required. The key 

On , step is to incorporate other readily available information about the system through the framework of a 

^^ • so-called generalized model. Our new approach may help to anticipate future catastrophic regime shifts 

in nature and society, allowing humankind to avert or to mitigate the consequences of the impending 
change. 
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Introduction 



Critical transitions are sudden, long-term changes in complex systems that occur when a threshold is 
crossed [1]. Many systems are known to be at risk of such transitions, including systems in ecology [2], 
climate research [3] , economics [1] , sociology [5] and human physiology [B] . Examples of critical transitions 
d ' in ecology include shifts in food web composition in shallow lakes [7], degradation of coral reefs [5], 

degradation of managed rangelands [3] , and desertification [TU] . 

Warning signals for impending critical transitions are highly desirable, because it is often difficult 
to revert a system to the previous state once a critical transition has occurred [2l[Tl]. If an accurate 
mathematical model of the system is available then critical transitions can be predicted straight-forwardly, 
either by numerical simulation or by direct computation of the dynamical thresholds. For real world 
complex systems, however, sufficiently accurate models are in general not available, and predictions based 
on models of limited accuracy face substantial difficulties HH . Recent research has therefore focused on 
model- free approaches that extract warning signals from observed time series [13] . Two of the most widely 
used approaches are increasing variance |14| and autocorrelation |15| , both of which are caused by critical 
slowing down J16j . Other approaches consider warning signals based on skcwness |17j, flickering |18j and 
spatial correlation }19j . 

In practice, robust calculation of early warning signals such as these is often limited by the need 
to acquire sufficient data. Especially in ecological applications, acquiring time series with the desired 



accuracy and sampling rate poses a significant challenge. One strategy for reducing the demand for time 
series is to utilize other information that may be available. This highlights the need for intermediate 
approaches, which can efficiently incorporate available information without requiring a fully specified 
mathematical model. 

In the present Letter, we propose an approach for the prediction of critical transitions based on the 
framework of generalized modeling |20|21| . The approach allows available information to be used to reduce 
the amount of time series data required or to increase the quality of the prediction. We demonstrate 
the applicability of the proposed approach by considering a simple one-population model, a previously 
studied fisheries model and a tri-trophic food chain. 

Methods 

Generalized modeling 

Suppose that a system has been identified as being at risk of a critical transition. Even if very little 
specific information is available, the dynamics can generally still be captured by a so-called generalized 
model |20| . Such a model captures the structure of the system, without restricting it to specific functional 
forms. 

To formulate a generalized model we first identify important system variables (say, abundance or 
biomass of the populations in the system) and processes (for example, birth, death, or predation). As a 
first step, the generalized model can then be sketched in graphical form, such as in Fig. [1] below. This 
graphical representation is sometimes called a causal loop diagram [22] . 

To obtain a mathematical representation of the model we write a dynamical equation for each variable 
Xi. In these equations we represent the processes by general functions. For instance we can model a 
single population Xi subject to gains G and losses L by an ordinary differential equation 

^Xi = G(Xi)-L(Xi), 

or as a discrete-time map 

X^,t+i = G{Xi,t)^L{Xi,t). 

Note, that in contrast to conventional models, we do not attempt to describe the processes G and L by 
specific functional forms. Instead, we use unspecified functions GO and L() as formal placeholders for 
the (unknown) relationships realized in the real system. 

Calculation of early warning signal 

We assume that before the critical transition, the system can be considered in equilibrium. We emphasize 
that this does not require the system to be completely static or closed in a thermodynamic sense, but 
that, on the chosen macroscopic level of description, the system can be considered at rest. Additionally, 
the system is subject to a slowly changing external parameter that puts it at risk of undergoing a critical 
transition. The system is therefore at equilibrium only on a certain timescale. In the following we refer 
to this timescale as the fast timescale, while the dynamics of the whole system, including the slow change 
of the external parameter, takes place on the slow timescale. 

Using the definitions above critical transitions can be linked to instabilities (bifurcations) of the fast 
subsystem [23] . For detecting these instabilities we construct the Jacobian matrix, a local linearization of 
the system around the steady state [23]. A system of ordinary differential equations (ODEs) is dynamically 
stable if all eigenvalues of the Jacobian have negative real parts, whereas a discrete time map is stable 
if all eigenvalues have an absolute value less than one. Critical transitions are thus signified by a change 



in the external parameter causing at least one of the eigenvalues to cross the imaginary axis (ODE) or a 
unit circle around the origin (map). 

To warn of impending critical transitions we monitor the eigenvalues of the Jacobian of the fast 
subsystem, which usually change slowly in time. A warning is raised if at least one of the eigenvalues 
shows a clear trend toward the stability boundary (for ODEs, zero real part; for maps, absolute value of 
one). The Jacobian itself can be computed directly from the generalized model, but will contain unknown 
terms reflecting our ignorance of the precise functional forms in the model. Previous publications j20j have 
shown that these unknowns can be treated as well-defined parameters with clear ecological interpretations. 
In the present applications we estimate the unknowns in the Jacobian matrix from short segments of time 
series data. Thereby, a pseudo-continuous monitoring of the eigenvalues of the fast subsystem is possible. 

We note that with given time series data estimating the generalized model parameters is simpler than 
estimating the entries of the Jacobian matrix directly, because the generalized model already incorporates 
structural information about the system. Further, many of the parameters in the generalized model may 
already be available in a given application, because they refer to well-studied properties of the species, 
such as natural life expectancy or metabolic rate. 

Results 

We applied the proposed approach to three case studies, focusing on a generic population with an Allee 
effect, a fisheries example, and a tri-trophic food chain. 

Simulation with Allee effect 

Alice effects, that is, positive relationships between per-capita growth rate and population size, are 
postulated in many populations and have been conclusively demonstrated in some |25j . A population 
with an Allee effect can suddenly transition from a stable, non-zero population size to unconditional 
extinction [26]. We constructed the Jacobian of such a system from a suitable generalized model (see 
Appendix SI in Supporting Information). This Jacobian contains unknown parameters that depend 
on the specific system under consideration. An early warning signal is obtained by estimating these 
parameters from time series data of the population size and birth rate, and monitoring the eigenvalues 
of the Jacobian as external parameters change. 

For testing the early warning signal, described above, we simulated a simple model 



d^_ AX 

dt k'^ + X 



-X = T^^^-mX + am, (1) 



based on Yeakel et al. [17] • We included a small additive noise term cr^(i), with standard deviation a 
and autocorrelation {^{t)^{t')) = S{t — t'). Throughout the simulation we slowly changed the mortality 
rate m according to m = 7.5 + 0.2i. A critical transition occurred, causing subsequent extinction of the 
population (Fig. [2]). 

The challenge we address is predicting the critical transition from a limited number (here, fifteen) 
of observations of population size and birth rate. We emphasize that we do not utilize any information 
on the functional forms of processes employed in the simulation, so that the prediction is based solely 
on the 15 observations and the assumed structural information (that is, one population subject to gains 
and losses). By estimating the parameters of the generalized model we determined the eigenvalues of the 
Jacobian as a function of time (Fig. Wp)- A clear increase in the eigenvalue is detectable well before the 
critical transition, giving ample warning of the impending collapse. 

Due to a phenomenon called bifurcation delay |23j , the population size did not start to change rapidly 
until well after (t w 13.5) the bifurcation point (f = 12.5). As previously observed by Biggs et al. [5S], 
management action to reverse the change in bifurcation parameter may successfully avert the critical 



transition even after the fast subsystem's bifurcation has occurred, if still within the range of the bifurca- 
tion delay. In the case of Fig. [IJb), the eigenvalue trend is directed more towards the last possible time 
that successful management action can be taken than towards the time of the actual bifurcation. 

Fishery simulation of Biggs et al. 

Our second case study focuses on an example from fisheries. Increased harvesting of piscivores can 
induce a shift from the high-piscivore low-planktivore regime to a low-piscivore high-planktivore regime 
[29] . Many fisheries are suspected to have undergone such transitions [30l|3T]. Based on the causal 
loop diagram (Fig. [1]), we formulated a discrete-time generalized model, describing the piscivore and 
planktivore populations at the end of each year, in the spirit of stock-assessment modeling (see Appendix 
SI). Thereby detailed modeling of the intra-annual dynamics was avoided. 

To test the warning signal, we generated time series data with a detailed fishery model by Biggs et 
al. [28], which was developed from a series of whole- lake experiments [32]. We note that this model 
differs significantly from our generalized model by a) accounting for the intra-annual dynamics and b) 
containing an additional state variable denoting the juvenile piscivore population. These discrepancies 
were intentionally included to reflect the limited information that would be available for the formulation of 
the generalized model in practice. In simulations the detailed model showed a transition to a low-piscivore 
high-planktivore regime as the harvesting rate was increased (Fig. [3^). 

From this simulation, we recorded the simulated piscivore and planktivore density and piscivore 
catch at the end of each year. Because the simulated data was very noisy we estimated the Jacobian's 
eigenvalues after smoothing the recorded data (see Appendix SI). In addition to the time series data, 
the information on the natural adult piscivore mortality and reproduction rate and the planktivore influx 
from refugia were required (see Appendix SI). This type of information can be reasonably well estimated 
for most systems. We confirmed that our predictions (reported below) are not sensitive to the specific 
values used. Indeed, a simple approach for estimating these parameters is to recognize that the initial 
state, before the critical transition, is stable. In a number of test trials we confirmed that any reasonable 
combination of parameters used that corresponded to an initially stable steady state provided an early 
warning signal comparable to the results reported below. 

An estimate of the Jacobian eigenvalues for the fisheries example is shown in Fig. [3j). As the system 
approaches the critical transition we observe that an eigenvalue approaches one, which signifies a critical 
transition for discrete time systems. This result is compared to the variance- based early warning signal 
computed by Biggs et al. [55] , which uses a much more densely sampled time series including intra-annual 
dynamics. The comparison shows that the approach proposed here produces a signal of similar quality 
(although possibly too early), while requiring significantly less time series data. Further, comparison with 
a variance-based signal using the same amount of time-series data as the generalized model shows that 
the generalized model-based signal is a much clearer early warning signal in this case. In particular, the 
variance signal only rises during or after the transition. 

Tri-trophic food chain 

For our final example we consider a tri-trophic food chain. In ecological theory food chains play a role both 
as a prominent motif appearing in complex food webs and as coarse-grained models. Using generalized 
models, a general Jacobian for a continuous-time model of the tri-trophic food chain can be derived (see 
Appendix SI and Gross et al. [33]). 

For generating example time series data we used a conventional equation-based model in which a 



producer biomass, Xi, predator biomass, X2, and top predator bioniass, X3, follow 

^Xi = A„Xi(if„ - Xi) - ^^ + ^i£.i{t) (2a) 

In these equations we assumed Holling typc-3 predator-prey interaction, logistic growth of the producer 
and linear mortality of the top predator and ^i{t) are additive noise terms with {£,i{t)£^j{t')) = SijS{t — t'). 
If any biomass decreased to zero, we suppressed the noise term so that the corresponding population 
remained extinct. 

We simulated the Eqs. ^ while increasing the mortality rate of the top predator. The resulting 
time series, for the chosen combination of parameters, show a slowly changing steady-state followed by a 
sudden transition to large oscillations, and a sudden cohapse of all three populations (Fig. |4|). 

To provide an early warning signal for the transition we recorded time series of the three biomasses 
and the top-predator's death rate, and estimated the parameters of the generalized model from smoothed 
segments of these time series. Even for the smoothed data we find that one of the eigenvalues is very noisy 
and sometimes positive. We believe that the presence of this eigenvalue reflects the response of the prey 
to fluctuations on the higher trophic levels and therefore exclude this value from our analysis. As m is 
increased toward the critical transition at the onset of oscillations, two eigenvalues show a clear increase 
toward zero real part (Fig.|4]). The two eigenvalues approach zero as a complex conjugate eigenvalue pair, 
which is indicative of the system undergoing a Hopf bifurcation |24j . This bifurcation generally implies 
a transition from stationary to oscillatory dynamics. The large oscillations combined with stochastic 
fluctuations then lead rapidly to extinction. The early warning signal constructed here thus indicates the 
transition to an oscillatory state signiflcantly before the transition occurs. 

Discussion 

In this Letter we proposed an approach for anticipating critical transitions before they occur. In particular 
we showed that generalized modeling of the system can facilitate the incorporation of the structural 
information that is in general available. 

We demonstrated the proposed approach in a series of three case studies. The first example showed 
that in simple systems even very few time points can be sufficient for clean prediction of the critical 
transition. The second example posed a hard challenge, where test data was generated by a model 
that differed considerably from the generalized model. Yet even in this case the generalized model 
significantly reduced the amount of data needed for predicting the transition. The third and final example 
demonstrated the ability of the proposed approach to distinguish between different types of critical 
transitions. 

In all case studies we found that the proposed approach can robustly warn of critical transitions in the 
presence of noise. We believe that the performance of the approach under noisy conditions can be further 
improved by subsequent refinements. Such refinements could include combination with dynamic linear 
modeling |14| . utilization of a parameter transformation (to 'scale' and 'elasticity' parameters) previously 
proposed for generalized models |20| . or the use of optimized sampling procedures. 

An advantage of the proposed approach is that it is generally becomes more reliable closer to a critical 
transition, where rates of change of state variables and other observables are generally larger, leading to 
better sampling. In contrast, methods based on variance or autocorrelation become statistically more 



difficult to estimate as the time series becomes less stationary. In this sense the proposed approach 
provides a tool complementary to established variance-based methods. 

An important assumption in our present treatment was that the dynamics of the fast subsystem could, 
at least at some level of description, be considered as stationary. Let us emphasize that this is not a 
strong assumption because even systems that arc primarily non-stationary, such as the fisheries example, 
can be modeled as stationary if a suitable modeling framework is chosen. Furthermore, ongoing efforts 
aim at extending the framework of generalized modeling to non-stationary dynamics, which may lead to 
a further relaxation of that assumption in the future |34j . 

In summary, we used generalized models to efficiently incorporate available information about a 
system without requiring detailed knowledge about that system. This led to the construction of early 
warning signals for critical transitions that required significantly less time series data than variance-based 
approaches. Thereby the proposed approach can contribute to the warning of critical transitions from a 
realistic sampling effort. 
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Figure 1. Schematic of the fishery knowledge that was incorporated into the generaUzed 
model. 
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Figure 2. Early warning signal for a single population with Allee effect, (a) Population time 
series (in relative units) generated by Eq. ([T]) (solid line, left axis), and yearly births (circles, right axis). 
The vertical dashed line indicates the time of the bifurcation. Parameters were A = 20, fc = 1, 
m = 7.5 + 0.2t and a = 0.1. The simulation was started at < < to allow for the decay of transient 
responses, (b) Eigenvalues estimated from the sampled data indicated with markers in (a), using the 
procedure described in the text. The eigenvalue was always real. 
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Figure 3. Early warning signals for the fishery simulation of Biggs et al. |28j . (a) Results for 
the adult piscivorc (blue line, left axis), juvenile piscivore (red line, left axis) and planktivore (green 
line, left axis) populations of the model of Biggs et al. [28], with the harvesting rate they specified of 
1.5 + 0.005i per year. From these results, data were sampled from A and F, once per year, after 
maturation and mortality for that year had been computed (black markers/thick black line). The 
annual piscivorc catch data used in the early warning signal calculation are also shown (thin black line, 
right axis), (b) Early warning signals: the eigenvalues (points, left axis) estimated by the generalized 
modeling approach; the intra-annual variance of the planktivore population (solid line, right axis) as 
calculated by the method of Biggs et al. [28] ; and the planktivore variance using only year-end data 
(dot-dashed line, right axis). Only one of the eigenvalues is visible at this scale. The eigenvalues were 
always real. The year-end variance warning signal was calculated using a sliding window of the previous 
15 year-end planktivore populations and a quadratic detrending within that window, which yielded 
generally stationary data for those windows that were before the critical transition. 
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Figure 4. Early warning signal for a critical transition in a tri-trophic food chain, (a) Time 
series of Xi (blue circles, left axis), X2 (green crosses, left axis) and X3 (red triangles, left axis) 
generated by Eq. ([2]). Only the data subsequently used in the early warning analysis were plotted. 
Some of the data during the oscillations between i = 40 and 45 were outside the scale of this graph, 
with X3 exceeding 30. The estimates of top- predator mortality used to calculate the early warning 
signal are also shown (black dots, right axis). Parameters were An = 4:, B ~ A, Ap — 2, Kp = 5, 
Kn = 10, K3 = 2, (cTi, 0-2, (T3) = (5, 0.02, 0.01) and m(i) = 0.4 + 0.02t. (b,c) Real and imaginary parts of 
eigenvalues estimated from the data in (a), using the procedure described in the text. Eigenvalues arc 
denoted by dots and circles; a dot within a circle indicates two eigenvalues had the same real value. The 
markers and colors used in (b,c) have no correspondence to those used in (a). A third (purely real) 
eigenvalue was not plotted, for reasons described in the main text. In (b), the horizontal dashed line 
indicates the stability boundary at zero real part, while the vertical dashed line indicates the time that 
the (Hopf) bifurcation occurred in the fast subsystem of Eqs. ([2|). 
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Appendix SI: Detailed method for early warning signal 
calculation 
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This Appendix details the generalized models, the assumptions behind the generalized models, and 
the numerical scheme used to generate the early warning signals for the three case studies in the main 
text. 

Single population with Allee effect 

In generating an early warning signal for critical transition in the population with Allee effect, we assumed: 
That the equation 



• 



• 



• 



• 



^^X = B(X)-M{X,fi) (3) 

is a good generalized model for the system, where B{X) and M{X,(i) are the birth and death 
rates of the population, respectively. Of these two processes, only the death rate is affected by 
the external parameter ^ pushing the system towards the critical transition. This factor could 
be anything that changes the death rate: the spread of a new disease, the appearance of a new 
predator, or habitat destruction, for example. 

Access to a time series of the population X{t). We refer to the population observations as AT,, taken 
at times i^, i = 1, . . . , n. 

Access to a time series of the birth rates Bi (the death rate Mi would also be acceptable). For 
simplicity, we assume these observations are also at times ti. 

That the mortality M{X,fi) is linear in X (although the coefficient of this linearity may change 
with fi). 



From Eq. ([3]) , the eigenvalue of the system's dynamics near its steady state is given by 

A = i3'(A)-M'(A,/i), (4) 

where the prime denotes the derivative with respect to X . This equation shows that the gradients B' and 
M' are sufficient to calculate A. Since the birth rate Bi and the population Xi have been directly observed, 
B'i = ABi/AXi could therefore be computed immediately, where we use the notation AGi = Gi — Gi-i. 
A discretization of Eq. (|3]) gives Mi = Bi — AXi/Ati. Because M depends on /i as well as X its derivative 
could not be estimated in the same way as B' . Using the assumption of linearity of M, we estimated 
M/ = Mi/Xi. (Suppose M = kX. Then M' = k = M/X.) We could then estimate the eigenvalue 

These derivative estimators are one-sided, in that they approximate /'(x) « [fix) — f{x — Ax)] /Ax. 
This one-sided estimation involved a loss in accuracy but allows the eigenvalues to be estimated at the 
most recent observation time, which is important when attempting to predict an imminent transition. 



12 

Fishery simulation of Biggs et al. 

We assumed the following knowledge: 

• The populations of adult piscivores Ai and planktivores Fi at the start of year i 

• The number of adult piscivores harvested Qi = Q{Ai) during year i 

• That adult piscivores die from other causes at a per-capita rate m 

• That adult piscivores predate on the planktivores 

• That there is a net movement of planktivores into (out of) the foraging arenas from (to) refugia |35| 
of Ri in year i 

• That per-capita, and in the absence of interactions with other modeled populations, r piscivores 
would be born and mature into adult piscivores the following year 

• That however juvenile piscivores are predated on by planktivores |29| and are also controlled, for 
example by accidental predation, by the adult piscivores |32| . 

This knowledge is represented schematically in Fig. 1 of the main text. 

We constructed a generalized model of this fishery by modeling only the year-end piscivore populations 
Ai and planktivore populations Fi with a discrete map. As a consequence we could avoid explicitly 
modeling the complications of intra-annual dynamics such as piscivore reproduction and maturation. We 
wrote: 

A,+i = (1 - m){A, - Q,) + rA, - C, (5a) 

F,+i = F, + R, - D,. (5b) 

The first term of Eq. ([Sa]) gives the number of surviving adult piscivores after harvesting {Qi) and 
other causes of mortality (m). The second term rAi models reproduction without losses and the third 
Ci = C{Ai,Fi) gives the predatory or controlling effects of planktivores and adult piscivores on the 
number of new adult piscivores. Since C depends twice on the adult piscivore population, first through 
the number of juvenile piscivores born and again through the controlling effect of the adults, we assume 
C is quadratic in Ai. 

In Eq. (j5bp . Ri models the net influx of planktivores into the foraging arenas. Since the planktivore 
population was initially small, we assumed there was negligible outflow of planktivores, and that the 
inflow depended on a refuge size that was constant. Therefore Ri = R was constant. [In the other 
scenario considered by Biggs et al. [28], where habitat destruction causes the refuge size to decrease, 
Ri would not be constant.] The remaining term Di = D{Ai,Fi) models all other factors changing the 
planktivore population, most importantly predation by adult piscivores. Since the planktivore population 
was initially small, we assumed D to be linear in Fi. 

We generated data to test this early warning signal with the previously published fishery model of 
Biggs et al. [28]. Their model is a hybrid discrete-continuous system and explicitly models the adult 
piscivore, juvenile piscivore and planktivore populations. The intra-annual dynamics are continuous, 
modeling: the harvest of adult piscivores; the predation and control of planktivores and juvenile piscivores, 
respectively, by adult piscivores; predation on juvenile piscivores by planktivores; and the movement of 
both planktivores and juvenile piscivores between the foraging arenas and refugia. At the end of each 
year discrete update rules were applied that controlled the mortality of adult piscivores, maturation of 
the surviving juvenile piscivores into adult piscivores, and birth of more juvenile piscivores. There was an 
additive noise term on the planktivore population dynamics. For further details of the model including 
its mathematical formulation see Biggs et al. j28| . 
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From the output of this simulation we recorded annual adult piscivore density, planktivore density, 
and piscivore catch. These observations were sufficient to calculate the Jacobian of Eqs. ([S]), 



J 






where we use the notation dH/dx = H'^'^\ The eigenvalues A of the Jacobian, which may be complex 
numbers, arc the solutions of det( J — XI) = where / is the (2 x 2) identity matrix. 

Our procedure was as follows. From the assumption of linearity in A, Q^ ~ Qi/Ai can be calculated 
immediately. The terms Ci and Di can be found by solving Eqs. ([5|). From the assumptions that C is 
quadratic in Ai and D linear in Fi, it follows that C^ = 2Ci/Ai and D^ = Di/Fi. Finally, by solving 
linear Taylor expansions of the form 

AC, = C,;(^) A^, + Cf^ AF, , (6) 

the derivatives C/ and D^ were estimated. 

This generalized modeling approach required estimates of the adult piscivore mortality m and repro- 
duction rate r, and planktivore influx R from refugia. We assumed the fishery manager could estimate 
m = 0.4 and r = 1.1 from knowledge of the fish species and the fishery. Effective values m = 0.5 and 
r = 1 can be derived from the actual simulation's equations |28j . The planktivore influx R would be 
harder to estimate in practice. In the absence of such knowledge, the ecosystem manager could make use 
of the fact that the fish populations were initially stable. This knowledge constrains R to the range 0.1 
to 1.5, as values outside these ranges produce an eigenvalue with magnitude greater than one between 
times i = to 10. We used R — 0.9. Indeed, this knowledge could also be used to constrain m and r to 
the ranges 0.1 to 0.7 and 0.95 to 1.5, respectively. We found that any combination of m, r and R that 
gave an initially stable system robustly predicted the later critical transition. 

Since the ecosystem approached the critical transition more slowly than the model with AUee effect 
above, to remove noise we found it necessary to smooth estimates of the gradients. For each new obser- 
vation at time ii, we fitted a second-order polynomial (since the method eventually results in differences 
of up to second order) to the record of each observed time series up to that time. Since it is the gra- 
dients at the current time that are most important, in the fit we applied a weight to each data point 
j of exp[(tj — ti)/T\ with T = 10. Gradients were then extracted directly from the coefficients of the 
polynomial fit. This simulated the calculation of early warning signals 'on the fly', as new observations 
become available. 

Tri-trophic food chain 

To obtain an early warning signal for this transition, we constructed a generalized model of the tri-trophic 
food chain 

^X^=A{Xi)-G{X^,X2) (7a) 

at 

^X2^G{Xi,X2)-H{X2,X3) (7b) 

at 

^Xs = H{X2, Xs) - M(X3, fi), (7c) 

at 

where Xi, X2 and X3 are the biomasses of the three species in the food chain, A{Xi) is the primary 
production rate of Xi, G{Xi,X2) and H{X2,X3) are the rates at which X2 and X3 consume Xi and 
X2, respectively, and M^X^, jj) is the total mortality per unit time of the top-level predator X^, with the 
argument /i denoting it is subject to a changing external influence. 
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We assumed: aceess to observations of all three biomasses Xi, X2 and X^, and the top-predator 
mortality M^X^^fi); that this mortality is linear in X3 so that M = m{^)X3; and that the predation 
rates G(Xi, X2) and H{X2, X3) are linear in their respective predator biomasses. These assumptions and 
observations were sufficient to calculate the Jacobian of Eqs. 0, and therefore the Jacobian's eigenvalues, 
by the same techniques as in the previous examples. The time series were first smoothed by the same 
filter as in the fishery simulation, but with r = 20. The approach could also easily be extended to the 
case of incomplete biomass conversion efficiencies, if estimates of those efficiencies were available. 
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